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ABSTRACT 

We have observed 6 clouds along the line-of-sight toward W49A using the 
Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based 
observatories. The ortho-H20 lio loi and OH (1665 and 1667 MHz) 
transitions are observed in absorption, whereas the low- J CO, ^^CO, and C^^O 
lines, as well as the [C I] ^Pi — ^ Pq transition, are seen in emission. The emission 
lines allow us to determine the gas density (n ~ 1500 — 3000 cm~'^) and CO 
column densities {N{CO) ~ 7.9 x 10^^ — 2.8 x 10^''' cm~^) using a standard 
Large Velocity Gradient analysis. 

By using both the o-Hg^O and 0-H2O absorption lines, we are able to 
constrain the column-averaged 0-H2O abundances in each line-of-sight cloud to 
within about an order of magnitude. Assuming the standard N(H2)/N(C0) 
ratio of 10^ we find N(o-H20)/N(H2) = 8.1 x 10"^ - 4 x 10"^ for three clouds 
with optically thin water lines. In three additional clouds, the H2O lines are 
saturated so we have used observations of the H2^0 ground-state transition to 
find upper limits to the water abundance of 8.2 x 10^^ — 1.5 x 10^^. We measure 
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the OH abundance from the average of the 1665 and 1667 MHz observations and 
find N(0H)/N(H2) = 2.3 x 10"^ - 1.1 x 10"^ The 0-H2O and OH abundances 
are similar to those determined for hne-of-sight water absorption features 
towards W51 and Sgr B2 but are higher than those seen from water emission 
hues in molecular clouds. However, the clouds towards W49 have lower ratios 
of OH relative to H2O column densities than are predicted by simple models 
which assume that dissociative recombination is the primary formation pathway 
for OH and H2O. Building on the work of Neufeld et al. (2002), we present 
photo-chemistry models including additional chemical effects, which can also 
explain the observed OH and H2O column densities as well as the observed 
H2O/CO abundance ratios. 

Subject headings: ISM:abundances ISM:individual(W49) ISM:molecules 
ISM:clouds starsrformation 



1. Introduction 

Since its launch in 1998 December, the Submillimeter Wave Astronomy Satellite 
(SWAS) has observed emission from the ground-state (lio Iqi) transition of ortho-H20 
(hereafter 0-H2O) in a number of molecular clouds (e.g. Melnick et al. 2000a; Snell et 
al. 2000a, 2000b). While it has been possible to determine the H2O abundances from the 
SWAS observations, these analyses depend on the often poorly contrained gas densities and 
temperatures. Observations of absorption from the 0-H2O ground-state transition, however, 
have the advantage that the cohimn density is simply proportional to the optical depth in 
the line. Therefore, the H2O abundance can be derived without precise knowledge of the 
temperature or density of the absorbing gas. 

In this paper we present observations of water absorption toward the molecular cloud 
complex W49A. W49 is 11.4 kpc from the Sun and 8.1 kpc from the Galactic Center 
(Gwinn et al. 1992). It consists of a supernova remnant (W49B) and an HII region 
(W49A) separated by ~ 12'. W49A has long been known as a site of extremely active 
star formation, due to its association with one of the most powerful H2O masers in the 
Galaxy (Genzel et al. 1978), and the fact that it is one of the most luminous regions in the 
Galaxy {L^oi ~ 10'^ L©; Ward- Thompson & Robson 1990). Physically, W49A is composed 
of a number of optically obscured, compact HII regions surrounded by a massive molecular 
cloud (M > 10^ M©; Mufson & Liszt 1977). W49A is broken up into three main IR peaks; 
W49SW, W49SE and, the strongest, W49NW (more commonly named W49N; Harvey et 
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al. 1977). The continuum emission from W49A has also been observed at 350, 800, and 
1100 //m by Ward-Thompson and Robson (1990) and at 1300 //m by Schloerb et al. (1987). 

Spectrally, W49A is extremely complex, containing numerous features contributed by 
W49A itself, as well as additional clouds associated with the Sagittarius spiral arm (which 
crosses the line-of-sight twice) . W49A and the associated hne-of-sight clouds have been the 
subject of numerous spectroscopic studies (both emission and absorption) involving species 
such as HI (Radhakrishnan et al. 1972), HCO+ and HCN (Nyman 1983), OH (Cohen and 
Willson 1981), H2CO (Mufson and Liszt 1977; Bieging et al. 1982), CO (Mufson and Liszt 
1977), SiO (Lucas and Liszt 2000), CS (Greaves and Wilhams 1994), and 01 (Vastel et al. 

2000) . 

The strong continuum provided by this source, and the well-studied deep absorption 
features in the gas associated with line-of-sight clouds, provide an excellent opportunity to 
measure the H2O abundances in a number of different molecular clouds. Similar analyses 
of the water abundances have recently been carried out for line-of sight clouds toward 
Sagittarius B2 ( Neufeld et al. 2000; Cernicharo et al. 1997), Sagittarius A* (Moneti et al. 

2001) , and W51 (Neufeld et al. 2002). 



2. Observations 

Since SWAS has two independent receivers, we were able to simultaneously observe 
the [C I] ^Pi ^3 (i, ^ 492.1607 GHz) and O2 3, 1 ^ 3, 2 (i/ = 487.249 GHz) transitions 
in the upper and lower sidebands of receiver 1, and the ^^CO J = 5 — > 4 (i/ = 550.926 
GHz) and 0-H2O lio — > loi {^^ — 556.936 GHz) transitions in the upper and lower sidebands 
of receiver 2. However, when receiver 2 was tuned to the o-H2^0 lio loi transition 
(u = 547.676 GHz), SWAS loses the abihty to observe 0-H2O and ^^CO J = 5 ^ 4. The 
integration times and la rms noise levels for each line are listed in Table 1. The SWAS 
beamsize at 557 GHz is 3.3' x 4.5' while at 490 GHz it is 3.5' x 5.0'. The main beam 
efficiency {i]mb) for SWAS is 90%. For more information about the SWAS instrument, data 
acquisition, and data analysis, see Melnick et al. (2000b). 

The SWAS observations of W49 were obtained in four different observing periods; 
1999 April and May, 1999 September and October, 2000 March - May, and 2000 
October. The central position (located at W49N) is given as a(J2000) = 19^10"13.5^ 
5(J2000) = 09° : 06' : 29". The reference position was ~ 1.5° from the central position 
(a(J2000) = 19'*15'"19.8^ 5(J2000) = 08° : 15' : 22") and was chosen based on an absence 
of ^^GO J = 1 ^ emission. 
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In order to understand the H2O and OH absorption, it is important to know the 
spatial extent of the submilUmctcr continuum. To this end, wc have mapped a region in 
W49 comparable to the SWAS beam at 450 fim using the SHARC camera at the Caltech 
Submillimeter Observatory (CSO). These observations were taken in 1999 September. 
The flux levels were calibrated using Uranus as the calibration source. In addition to the 
SHARC maps, we also have data from the on-board SWAS continuum detectors in the 
fast-chop (4 Hz) observing mode (see Melnick et al. 2000b for details). 

Wc observed ^^CO J = 1 ^ (z/ = 115.271 GHz), ^^CO J = 1 ^ (z/ = 110.201 GHz), 
and C^^O J = 1 ^ {u = 109.782 GHz), at the 14m Five Colleges Radio Astronomy 
Observatory (FCRAO) using their 32 element array receiver (SEQUOIA). The ^^CO 
J = 1 — > observations were taken in 2001 August on an 8 beam x 8 beam grid with 
44" spacings (the resolution of the FCRAO antenna). The ^^CO and C^^O J = 1 ^ 
observations were obtained in 2002 November. For comparison with the SWAS H2O 
observations these data were smoothed to a resolution of 1.5' to match the size of the 
continuum source in W49 (see Section 3.1). The main beam efficiency {rj^b) of the FCRAO 
is 49% 

Observations of the ^^CO J = 3 — > 2 transition were made at the KOSMA (Kolner 
Observatorium fiir Submm-Astronomic) 3m telescope near Zermatt, Switzerland. The 
observations were taken on 2000 November 19 using the MRS (medium resolution acousto 
optical spectrometer) backend with a resolution of 360 kHz, corresponding to 0.31 km s~^ at 
345 GHz. The beamsize at 345 GHz is 82" and the main beam eflicicncy (r^mb) is 78%. For 
comparison with the SWAS observations, these data were smoothed to match the spatial 
resolution of the SWAS ^^CO J = 5 ^ 4 observations. 

In order to further examine the line-of-sight clouds towards W49A, observations of the 
F = 1 ^ 1 (1665 MHz) and F = 2 ^ 2 (1667 MHz) transitions of OH were obtained at 
the Arccibo Observatory in 2001 December. At these frequencies, the Arecibo beam size 
is approximately 2.6' x 3' (FWHM) and has a main beam efficiency at 1665 & 1667 MHz 
of 80%. We also obtained observations of the 1612 and 1720 MHz transitions of OH at 
the Arecibo observatory in 2003 October allowing us to investigate LTE versus non-LTE 
excitation in the OH lines. 

The on-source integration times and 1 a rms noise levels for each of the observed 
spectral lines are presented in Table 1. All spectra in this paper are presented in units of 
T^. For subsequent analysis we use the main beam temperature (T^^;,) which is defined by 
the standard equation T^b — TX/rjmb where rjmb is the main beam efficiency. 
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3. Results 

Figure 1 shows the 450 fim continuum emission obtained using SHARC at the CSO. 
The continuum emission is strongly concentrated on W49A itself but with weaker tendrils of 
emission extended approximately 0.75' from the center. We use the 450 /xm map of W49 to 
confirm the continuum level in the SWAS 0-H2O and o-H2^0 spectra. First, we extrapolate 
the integrated 450 iim continuum emission to 540 iim (the wavelength of the SWAS 0-H2O 
line) assuming a power law of the form S oc A~^'^. Then, we convert the integrated fiux at 
540 /xm to an equivalent main beam brightness temperature (Tmb) through the standard 
equation: S = ^^"^'g''^ . Using parameters appropriate for the SWAS beam at 540 fim, we 
obtain T„ib ~ S{Jy) /13982. The integrated fluxes for W49 are 1827 Jy at 450 /im and 965 
Jy at 540 /xm, corresponding to Tmb = 0.07 K at 540 /xm. 

At 539 fim, the SWAS on-board continuum detectors yield a Tmb of 0.1 K for W49A. 
The fact that the SWAS results are ~ 30% higher than SHARC observations probably 
reflects uncertainties in the extrapolation to 540 fira, atmospheric effects, and the fact that 
the SHARC maps do not perfectly match the SWAS beam size, shape, and orientation. 
It should be noted, however, that the SWAS measurements may not be entirely accurate, 
due to a relatively slow chop rate. Therefore, for the subsequent analysis in this paper we 
instead use the fact that four of the H2O absorption features are at the same depth to within 
the noise (the ~ 18 km s~^ component from W49 itself, and the 39.5, 59.6, and 63.3 km 
features; see Figure 2). This strongly suggests that these fines are saturated and thus, the 
continuum flux level can be set by the depth of these four features. The average depth of 
these four lines is 0.08 K (T|) which agrees satisfactorily with both our measurement from 
SWAS' on-board continuum detectors and our extrapolation of the SHARC data. 

Figure 2 shows the 0-H2O lio loi, [C I] ^Pi -^^ Pq, and ^^CO J = 5 4 observations 
from SWAS. Velocities less than ~ 30 km s~^ correspond to emission/absorption intrinsic to 
W49A and these features are not further analysed in this paper. At velocities greater than 
~ 30 km s~^, Figure 2 shows six distinct 0-H2O absorption dips at LSR velocities of 33.5, 
39.5, 53.5, 59.6, 63.3 and 68 km s"^ Gaussian fits to the baseline subtracted spectra are 
fisted in Table 2 for each of the water absorption features. The same six features are seen 
in emission in the [C I] spectrum (Figure 2), although at sfightly different LSR velocities. 
The sfight offset between the [C I] and 0-H2O fine centers may be partly due to the high 
optical depths of the water lines. The Gaussian fits to the baseline subtracted carbon fines 
are given in Table 3. No ^'^CO J = 5 — 4 lines were detected at velocities > 30 km s"-*^. 
The O2 transition was not detected to a 3a level of 15 mK. 

Figure 3 shows the ^^CO, ^^CO and C^^O J = 1 ^ observations from FCRAO and 
the ^^CO J = 3 — > 2 observations from KOSMA. All spectra have been convolved to a 
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circular, Gaussian beam of 3.8' to match the beam resolution of the SWAS ^^CO J = 5 — > 4 
observations . Figure 3 shows that while all six H2O absorption features are seen in both 
CO lines and the ^^CO J = 1 — lines, the less abundant C^^O J = 1 — > only shows 
the 3 strongest components. It should be noted that, while the spectra in Figure 3 have 
been smoothed to match the spatial resolution of the SWAS beam, for our subsequent LVG 
analysis (Section 3.1), we use data smoothed to the size of the background continuum 
source (1.5'). Since we are trying to estimate the CO column density in the gas which is 
absorbing the H2O, it is more relevant to coadd only those ^^CO and C^^O spectra that fall 
in front of the continuum source. Therefore, we coadd the ^^CO and C^^O spectra within 
~ 0.75' of the central position (Figure 1). Gaussian fits to the detected spectral features in 
the data smoothed to 1.5' (with first order basehnes subtracted) are given in Table 4. 

Figure 4 plots the hne to continuum flux ratio for the SWAS observations 
of both 0-H2O and o-H2^0. The hne to continuum flux ratio is given by 

Fi/Fc = [ATynbiline) + Tm,b{cont)]/Tmb{cont) where T„ib{cont) is the measured main 
beam, SSB continuum flux level (see above) and AT„ih{line) is the baseline subtracted 
antenna temperature (ATJ^; listed in Table 2) divided by the main beam efficiency. 
Although there are no clear o-H2^0 lines at velocities > 30 km s^^, we provide upper limits 
in Table 2 by fitting Gaussians with LSR velocities and FWHM linewidths fixed at the 
values provided by the [C I] emission lines. We used the [C I] spectrum as a template, rather 
than the 0-H2O spectrum, due to the large optical depths of the water hne. Therefore, 
Table 2 provides hmits to the o-E^O hne strengths and integrated intensities. Table 2 also 
lists the line to continuum flux ratios (seen in Figure 4) and the line center optical depth 
{To^-ln{Fi/F,)). 

Since the continuum level is set by the average depth of the four deepest absorption 
features, the 39.5, 59.6, and 63.3 km s~^ 0-H2O features in Figure 4 reach line-to-continuum 
values close to, or less than, zero (and are, in fact, as strong as the absorption feature 
associated with W49A itself). Therefore, these features are most likely saturated, and 
we would consequently derive a null or negative line flux and an unphysical opacity. For 
these three lines, we provide a minimum line center opacity in Table 2 by adding the la 
error of the Gaussian fits ( T\ = 0.006 K) to the measured continuum level. The other 
three features (at 33.5, 53.5, and 68 km s^^) are weaker and appear optically thin. It is 
important to note however, that while we assume these components to be optically thin 
they may, in fact, have a higher optical depth. The relative hne strengths of the six ^^CO 
J = 1 — > components show significant variation over the region mapped, implying the 
presence of cloud structure on scales smaller than the SWAS beam. Therefore, we cannot 
rule out the possibility that the 33.5, 53.5, and 68 km s~^ features are actually saturated 
and only appear to be optically thin because they do not completely cover the continuum 
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source. Nevertheless, for the subsequent analysis in this paper we will assume that these 
three components are unsaturated. 

Figure 5 plots the hne to continuum flux ratio for the Arecibo observations of the 
OH 1665 MHz and 1667 MHz transitions. Since the front end amplifiers at the Arecibo 
observatory have uniquely defined pass bands, there is no single sideband to double sideband 
conversion issue. Thus, the measured SSB continuum levels {Tmb{cont)) are 311.3 K at 1665 
MHz and 315.3 K at 1667 MHz. All six of the features seen in emission in [CI] are seen in 
absorption in both OH lines. Tabic 5 lists the Gaussian fit parameters to all six components 
(with a first order baseline removed), as well as the line to continuum flux ratios and 
associated line center optical depth. The lines at 39.3, 54.2, 62.9, and 68.8 km s^^ are all 
within 25% of the 5/9 integrated intensity ratio expected for the 1665/1667 lines in local 
thermodynamic equilibrium (LTE). The 1665/1667 ratios in the 33.6 and 59.8 km s~^ lines, 
however, are ~ 1/3 which is ~ 40% lower than expected for LTE. The weakness of these 
lines suggests that they are optically thin. Therefore, the 1665/1667 ratio of 1/3 seen in 
the 33.6 and 59.8 km s~^ lines is probably a result of non-LTE excitation and not a line 
saturation effect. In addition, our most recent OH observations (at 1612, 1665, 1667, and 
1720 MHz) taken at the Arecibo observatory strongly suggest that the 39.3 and 62.9 km 
features are also affected by non-LTE excitation. Therefore, we will restrict further analysis 
of the OH absorption lines to the 54.2 and 68.8 km s~^ components. 



3.1. CO and C" Abundances 

To determine the line-of-sight average H2O and OH abundances, we flrst need to 
estimate the CO column density corresponding to each of the six water absorption features. 
To do this, we use a Large Velocity Gradient (LVG) code to simultaneously fit the ^^CO 
J = 1 ^ 0, 5 ^ 4, and C^^O J = 1 ^ observations. Since ^^CO J = 5 ^ 4 emission 
was not detected in any of the clouds, we use the 1 a noise level in the spectrum as an 
upper limit to the line strength in order to better constrain the fitting routine. We do not 
include the ^^CO observations in this analysis due to the higher opacity in these transitions 
and chose to focus instead on the optically thin tracers. An investigation into the effect 
of ignoring the ^^CO data in the LVG fits shows that the derived densities and column 
densities change by less than a factor of 2. We also use the data listed in Table 4, which 
have been smoothed to the size of the background continuum source rather than to the 
spatial resolution of the SWAS beam. The densities and column densities determined from 
this data set, however, do not differ significantly from an identical analysis using data 
coadded over an entire SWAS beam. 
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We created a 20 x 20 grid of models in density-column density parameter space 
using a constant kinetic temperature of 8 K (Vastel et al. 2000), consistent with what 
is derived from our ^^CO observations assuming that the emission is optically thick and 
thermalized. The densities range from 1000 — 5000 cm~^ and the ^'^CO column densities 
(per velocity interval) range from 10^'^ — 10^^ cm~^/km s~^. To fit the observations of the 
^^CO and C^^O isotopologues we assume isotopic abundance ratios (with respect to CO) 
of 55 and 500 respectively. The observed line intensities are fit to the grid of LVG models 
using a minimization routine to find the density and column density combination that 
best fit the observations. Results from the LVG fitting are listed in Table 6 along with 
the estimated visual extinctions where we have assumed that N(H2) ~ 8 x VS^'^Ay and 
N(H2)/N(C0) ~ 10^. It is important to note that this only estimates the dust extinction in 
the CO emitting region and ignores any additional extinction which may arise in overlying 
layers of HI or H2. To check the sensitivity of our results to our assumed value of kinetic 
temperature, we performed an identical LVG analysis using kinetic temperature of 15 K 
instead of 8 K. The resultant densities change by less than a factor of three and the column 
densities by less than 25%. 

Using the H2 density in Table 6 we also calculate the neutral carbon (C°) column 
density from the SWAS [C I] — >^ Pq observations (using the same LVG code and again 
assuming a kinetic temperature of 8 K). The results are also tabulated in Table 6 along 
with the C°/CO abundance ratio. The C°/CO abundance ratio varies from 0.6 for clouds 
with inferred visual extinctions (Ay) of 1.9, to C°/CO = 3.8 for clouds with inferred Ay 
= 0.1. These numbers are consistent with C°/CO ratios seen at the low column density 
edges of UV illuminated giant molecular clouds (e.g. Plume et al 1999) and with previous 
observations of high-latitude and translucent clouds (Ingalls et al 1997; Stark et al 1996; 
Stark & van Dishoeck 1994). 



The SWAS ortho-H20 column densities are derived from a "curve of growth" analysis 
of SWAS o-H2^0 and 0-H2O spectra. The column density in the lower state (A^/)can be 
expressed as: 



3.2. H2O and OH Abundances 



gi Snu^ y/TT 1 

Qu c3 [1 - e-Wfc^-] 



AV{FWHM) 




(1) 
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where Tq is the optical depth at hne center. For the 0-H2O and o-H2^0 absorption hnes in 
W49, Tex is probably much less than hu/k — 27 K, and so equation (1) reduces to: 

_ ^Stt AVjFWHM) 

Using parameters appropriate for the SWAS observations: A(o — Hl^O) — 547.39//m, 
A(o - H2O) = 538.66/xm, A^^i = 3.5 x IQ-^ g-^ (for both 0-H2O and o-H^^O), and qi/qu = 1, 
the expressions for the column density reduce to: 

Ni{o - Hl^O) = 4.66 X 10^2 ro{H^^O)AV{FWHM), (3) 

and 

Ni{o - H2O) = 4.89 X 10^2 ro{H20)AV{FWHM), (4) 

where AV{FWHM) is in units of km s~^. Again, since the excitation temperatures in 
these clouds is low, A^; is, to first approximation, equal to the total column density. To is the 
line center optical depth of the particular molecular transition. 

Table 7 presents the ortho-water column density in each of the six line-of-sight 
absorption features towards W49A. The first row lists upper limits on TqAV for o-H2^0 
taken from the o-H2^0 observations listed in Table 2. The second row presents upper 
limits on the o-H2^0 column densities as calculated from equation (3). The third 
row gives the 0-H2O column density upper limit calculated under the assumption that 
N{H20) = 500 X N{H^^O). The fourth row lists ToAV for 0-H2O taken from the 0-H2O 
observations listed in Table 2. The fifth row presents the 0-H2O column densities as 
calculated from equation (4). The sixth row hsts the CO column density derived from our 
observations of CO and its isotopologues. The seventh row of Table 7 presents the 0-H2O 
abundance as determined from the o-H2^0 observations (i.e. by dividing row 3 by row 6) 
and assuming that the CO abundance (relative to H2) is 10~^. Row eight provides the 
same 0-H2O abundance except that in this case we use the 0-H2O observations directly. To 
obtain the total H2O abundances, the numbers listed in Table 7 must be multiplied by 4/3 
to account for the assumed ortho to para ratio. 

The o-H2^0 observations allow us to obtain upper limits for the 0-H2O abundance, 
whereas the 0-H2O observations set lower limits. Thus, by using both the o-Hg^O and 
0-H2O observations, we are able to constrain the 0-H2O abundances in each line-of-sight 
cloud to within about an order of magnitude. Table 7 shows that the water abundance in 
the three clouds with optically thin water hnes is between 8.1 x 10~^ and 4 x 10~^ . In 
the three additional clouds, where the H2O lines are saturated, the upper limits to the 
water abundance are between 8.2 x 10~^ and 1.5 x 10~^. These abundances are similar to 
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those determined for a line-of-sight ortho-water absorption feature towards W51 (Neufeld 
et al. 2002), Sgr B2 (Neufeld et al. 2000; Cernicharo ct al. 1997), and Sgr A* (Moneti et 
al. 2001). The water abundances in these absorption features are, however, higher than 
those seen from ortho-water emission hues in other giant molecular clouds (2 — 10 x 10~^; 
Snell et al. 2000b). An explanation for the low water abundances observed in typical giant 
molecular clouds observed by Snell et al. (2000b) is given by Bergin et al. (2000) who 
suggest that, in the dense interiors of these well-shielded cores, H2O can readily freeze out 
onto dust grains. 

To calculate the OH abundances, wc again use cquation(l). However, for OH, 
hiy < kTgx even in the cold absorbing clouds and, therefore: 

^_9l87n^^ AVjFWHM) 

hc^ A^i '^^^'^ 2VM ■ ^ ' 

Inserting the appropriate constants for the 1665 MHz (^liees = 7.11 x 10~^^ and 
gi/g^ = 3/3 = 1) and the 1667 MHz lines (Aieey = 7.71 x 10"^^ and gi/gu = 5/5 = 1) we 
obtain the following equations for column density: 

16 

Ntot{lG65) = — X Ar;(l665) = 4.30 x 10^^ Te^ To{1665)AV{FWHM), (6) 

and 

16 

Ntot{lGG7) = ^ X A^i(1667) = 2.38 x 10^^ Te^ To{1667)AV{FWHM), (7) 
5 

where again AV{FWHM) is in units of km s~^. The factors 16/3 and 16/5 account for the 
relative populations in the four OH hypcrfine levels. 

To determine the OH column densities, therefore, we need to know the excitation 
temperature. Since the lines at 54.2 and 68.8 km s^^ appear to be in LTE, we assign T^x 
= 8 K as we did for the CO analysis. The other four clouds, however, appear to suffer 
significantly from non-LTE excitation and, since we have no way of knowing the excitation 
temperature in these clouds, we cannot accurately calculate their column densities and do 
not consider them in the following discussion. Fortunately, the 54.2 and 68.8 km s~^ clouds 
are also those in which the H2O lines are optically thin, and so provide us with our best 
measure of the water column density. Table 8 shows that the OH column densities in the 
remaining two clouds (where N(OH) is the average of the 1665 and 1667 MHz results) range 
from 3.5 x 10^^ cm^^ to 1.4 x 10^^ cm^^, similar to the value of 8 x 10^^ cm~^ seen in the 
hue of sight cloud towards W51 (Neufeld et al. 2002). 
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4. Discussion 

4.1. H2O and OH Abundances: A New Perspective on Branching Ratios 

The H2O/OH abundance ratio can provide insight into the chemical networks that 
produce oxygen-bearing molecules. Recently, Neufeld et al. (2002) made predictions for 
the ratio of OH and H2O column densities by using a simple analytic model in which OH 

and H2O formation, through dissociative recombination of HaO"*" (producing either O, OH, 
or H2O), was balanced by photodissociation. Neufeld (2002) compared the results of the 
simple analytic model with the detailed results of a diffuse cloud PDR model (e.g. Kaufman 
et al. 1999) and found that the simple analytic model predictions were robust for clouds of 
gas density n = 100 cm~'^, for a range of cloud extinctions and cosmic ray ionization rates. 
Those model results show that cosmic ray ionization of H2 leads directly to the formation of 
H3O+ which is subsequently converted to H2O and OH through dissociative recombination. 

The OH and H2O abundances, however, depend sensitively upon the poorly constrained 
chemical branching ratios (/o, foH, and Jeqc) for dissociative recombination of HaO"*". The 
flowing afterglow laboratory experiment by Williams et al. (1996) found fn20'-fon'-fo — 
0.05:0.65:0.3. Two other experiments using a different technique, one using the ASTRID 
heavy-ion storage ring (Vejby-Christensen et al. 1997; Jensen et al. 2000) and the other 
using the CRYRING heavy-ion storage ring (Neau et al. 2000), yielded fn^o'-fou'fo = 
0.25±0.01 : 0.74±0.02 : 0.013±0.005 and fn^o-foH-fo = 0.18±0.05 : 0.78±0.08 : 0.03±0.06 
respectively. Observations of diffuse clouds in the ISM have yielded similarly discrepant 
results. Recent UV observations with HST towards the moderately reddened star HD 
154368 (Spaans et al. 1998) found a 3cr upper limit on /j^jO of 0.06, which is consistent 
with the flowing afterglow experiment. However, observations of a line of sight cloud 
seen in absorption against W51's continuum are, instead, consistent with the ASTRID 
and CRYRING experiments (/ffjO ~ 0.25; Neufeld et al. 2002). Our results, however, 
are consistent with neither of the experimental values for the branching ratios. In Figure 
6 we show the results of the simple analytic model for branching ratios /h20-/oh:/o of 
0.05:0.65:0.30 and 0.25:0.75:0.0 (the top and bottom dashed lines respectively). We also 
show the measured values (or upper/lower limits) for two of the six observed absorption 
features (for which we have OH column densities). As may be seen in Figure 6, the 
observations lie below and to the right of the line relevant for the branching ratio 
0.25:0.75:0.0, and are even less consistent with the smaller H2O branching ratio, at least in 
hght of the conclusions of the simple analytic model. 

Neufeld et al. (2002) provide an explanation for the discrepancy in the branching 
ratios. When formation of H2O and OH via dissociative recombination of HaO"*" is the only 
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formation route, the predicted ratio of column densities depends only on the branching 
ratios. If, however, there are significant other routes to OH or H2O, then the ratio of column 
densities will depart from the simple analytic model predictions. If even a small fraction 
of gas in the clouds is warm (above ~ 300 K) then neutral-neutral reactions, which have 
relatively high activation barriers, can begin to contribute to the production of OH and 
H2O. Neufeld et al. (2002) show that, by varying the temperature of the gas, it is possible 
to match either of the experimental branching ratios. This model, however, requires a small 
fraction of gas with temperatures in excess of 600 K to explain the OH and H2O column 
densities in the clouds observed towards W49 (Figure 6). In the remainder of this section, 
we present an additional explanation for the observed OH and H2O column densities which 
builds on the Neufeld et al. (2002) results. 

To investigate additional chemical effects on the column density ratio, we have run 
models similar to those in the Neufeld et al. (2002) paper, but extending to higher gas 
densities. In all of the models, we assume that a diffuse cloud is illuminated from both sides 
by a total interstellar field equal to Go=1.7 which corresponds to the current best estimate 
of the local interstellar radiation field (a value Go=l corresponds to the Habing (1968) UV 
field; 1.3 x lO^"' erg cm^^ s^^ sr^^). In accordance with the densities and column densities 
derived from our LVG analysis (Table 6), we compute the total water and hydroxyl column 
densities through clouds with gas densities n = 100, 10^ and 10^ cm~^ and with total visual 
extinctions of Ay — 1, 2, 3, and 4. The gas temperature is solved for self-consistently 
and, at the center of the clouds, is found to be 22 K for Ay = 2, 12 K for Ay = 3, and 
8 K for Ay — 4. These values of Ay are consistent with those estimated from our ^^CO 
observations since they are the total Ay through the cloud, whereas the values listed in 
Table 6 only consider the CO emitting region. In the surface layers (to Ay ~ 1 on each side 
of the cloud) there is essentially no CO (the CO/H2 abundance ratio < 10^''^) but the dust 
in these H and H2 layers still contributes to the overall extinction. Therefore, one needs to 
add approximately 2 magnitudes of extinction to the values listed in Table 6 to compare 
with our models. 

The destruction of OH and H2O need not be dominated by photodestruction as in the 
simple model of Neufeld et al (2002). OH can be destroyed by neutral- neutral reactions 
with atomic O. The photodestruction rate per OH molecule is proportional to Go, which 
is held fixed, but the destruction by O is proportional to n(0), the density of atomic O. 
Thus, at higher n, the latter mechanism dominates, OH is destroyed more rapidly and the 
H2O/OH column density ratio increases relative to the simple toy model. Similarly, at 
high extinction, the photodestruction rate decreases due to dust attenuation of the FUV 
field, but the neutral rate is unaffected, leading to a relatively higher rate of destruction 
of OH than H2O and an increased H2O/OH ratio. These effects can be seen in Figure 6 
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which shows an increase in the ratio of H2O to OH column density as the gas density and 
visual extinction increase. The fact that all four features seen in W49 lie at or below the 

= 0.25 curve (and much below the fH20 = 0.05 curve) suggests both that we can rule 
out a branching ratio of /i^jO = 0.05, and that a simple model which incorporates only 
dissociative recombination is inadequate. 

Based on ISO observations of the [O I] 63 fim. transition, Vastel et al. (2000) have 
suggested that CO is depleted by at least a factor of 6 in these clouds. If true, then our 
H2O and OH abundances would also need to be lowered by a similar amount. However, 
there are several potential difficulties with this interpretation. First, while the authors 
analysed their data as carefully as possible, it is inherently difficult to compare the 63 /xm 
[O I] absorption feature to the HI and molecular features due to the poor spectral resolution 
of ISO {/SVpwHM ~ 44 km s~^). Second, the HI column densities quoted by Vastel et al 
(2000) are a few x 10^^ cm~^ and the HI is fairly optically thick (r ~ 2 — 4; Lockhart & 
Goss 1978). Therefore, it is possible that they are underestimating the HI column density, 
in which case the atomic gas may account for more of the observed [O I] 63 //m absorption. 
Finally, if CO is freezing out on grains, then oxygen must suffer the same fate. Calculations 
by Bergin et al. (2000) show that atomic oxygen depletes onto grains even more readily 
than CO. In fact, in the line-of-sight clouds towards W49, our pure gas-phase models are 
able to produce results that are consistent with the observed H2O/CO column density 
ratios, without having to resort to freeze-out of CO molecules onto dust grains. In figure 
7 we present the results of our PDR calculations along with the observed H2O and CO 
column densities. The H2O/CO abundance ratio in the three clouds with optically thin 
water lines is between 8.1 x 10~^ and 4 x 10~^ and, in the three clouds where the H2O lines 
are saturated, the upper limits to the H2O/CO ratio is between 8.2 x 10^"^ and 1.5 x 10~^ 
(Table 7). Figure 7 shows that the models are consistent with the observations for Ay > 3 
and n > 100 cm^^. 

If depletion does play a role in these clouds, then the freeze-out of atomic oxygen may 
help to explain the observed H2O/OH abundance ratios in clouds with pressures more 
closely matched to those observed in other diffuse interstellar clouds. For example, the 
54 km s"^ feature is consistent with models having Ay between 3 and 4 and intermediate 
densities (n < 10"^ cm^"^). The 68 km s^^ feature, however, requires higher density models. 
The resultant pressures {nT ~ few x 10^ cm~^ K) in this cloud, therefore, is quite high for 
a diffuse interstellar cloud. However, it is possible to lower the gas density and still obtain 
high H2O/OH column density ratios if we include chemical reactions on the surfaces of 
grains. Preliminary calculations indicate that the inclusion of grain chemistry, involving 
O and C bearing species, leads naturally to the observed high H20-to-OH column density 
ratios. In these models, photodesorption of water ice formed on grain surfaces is the main 
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source of gas phase H2O throughout the clouds (atomic oxygen sticks to the grains and 
rapidly is converted to water ice which is photodesorbed) . This results in a higher water 
column density than found in the models without grain chemistry. OH, on the other hand 
continues to be formed primarily via the dissocative recombination of H3O"'". The point 
labeled as "grain" in Figure 6 shows the effect of adding grain surface chemistry even in 
very low density gas (n ~ lO^cm"^). Although the inclusion of grain surface chemistry is 
not needed to explain the H2O/OH ratios in the clouds towards W49, the attractive feature 
of this model is that we can still produce high H20-to-OH column density ratios in lower 
density gas (n < 10^ cm~^) with pressures more reflective of diffuse clouds. We will explore 
the effects of grain chemistry in diffuse clouds in more detail in a subsequent paper. 

4.2. CO & [CI] Intensities 

We have also used our models to predict the strengths of individual [C I] (492 GHz) 
and ^^CO J = 1 — > emission hues. Since [C I] is observed in emission with a relatively 
large beam, we need some way to estimate the emission produced by the same gas that 
produced the H2O and OH absorption features. To do this we compare the ^^CO intensity 
averaged over the ~ 1.5' continuum source (Table 4) to that averaged over the SWAS beam 
size; the ratio of these two intensities is then used to scale the observed [C I] emission 
in order to estimate how much of it comes from the direction of the continuum source. 
Corrected values for [C 1] are given in Table 3. 

In Figure 8, we plot the observed (corrected) [C I] and ^^CO integrated intensities 
and the intensities predicted from the standard FDR models. The uncorrected [C 1] and 
^^CO integrated intensities smoothed to the resolution of the SWAS [C 1] observations 
are not shown, but do not differ significantly from the points plotted in Figure 8. The 
observed intensities are not well matched by the models. For instance, the models can only 
match the observed CO intensities by resorting to high gas densities, a known problem 
with FDR models noted by other authors (e.g. Bensch et al. 2003). These high density 
models, however, overpredict the [C I] intensity. There are a variety of ways in which this 
shortcoming could be resolved. For instance, if CO self-shielded more effectively than is 
generally assumed, conversion of C° to CO would happen closer to the cloud surface. Then 
the CO intensity would be higher and the [C I] intensity would be lower. Another solution 
might be the use of constant pressure FDR models; our models assume that the gas density 
is constant. If the [C I] emission came from gas with a density below the critical density for 
the 492 GHz transition, while CO emission came from higher density (though cooler) gas, 
then the desired effect might also be achieved. However, in a preliminary constant pressure 
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calculation with low surface density and low UV field, the density contrast between the 
[C I] and CO regions was only a factor ~ 2; not enough to explain the observed differences. 
The discrepancies between the [C I] and ^^CO observed and model line strengths will be 
investigated more thoroughly in a future paper. 



5. Conclusions 

We have analysed emission and absorption lines from 6 clouds along the line-of-sight 

toward W49A. Using the Submillimeter Wave Astronomy Satellite (SWAS), we observed 
[C I] ^Pi ^3 Po (z^ = 492.1607 GHz), ^''^CO J = 5 ^ 4 (z/ = 550.926 GHz), 0-H2O 
lio ^ loi {r^ = 556.936 GHz) and o-H^^O lio ^ loi transition (z/ = 547.676 GHz). We 
observed ^^CO J = 1 ^ {u = 115.271 GHz), ^^CO J = 1 ^ {u = 110.201 GHz), and 
C^^O J = 1 ^ (z/ = 109.782 GHz), at the 14m Five College Radio Astronomy Observatory 
(FCRAO) and the ^^CO J = 3 ^ 2 transition at the KOSMA (Kolner Observatorium 
fiir Submm- Astronomic) . We also observed the 1665 and 1667 MHz transitions of OH 
at the Arecibo Observatory, and mapped the 450//m continuum emission at the Caltech 
Submillimeter Observatory (CSO). 

The 0-H2O and OH (1665 and 1667 MHz) transitions are observed in absorption, 
whereas the other lines are seen in emission (with the exception of o-H2^0 and ^^CO 
J = 5 ^ 4 which were not detected). Using the emission fines of ^^CO and C^^O we derive 
gas densities of 1500 — 3000 cm~^ and CO column densities of ~ 7.9 x 10^^ — 2.8 x 10^'^ 
cm~^ via a standard Large Velocity Gradient analysis. The observations of 0-H2O and OH 
absorption have the advantage that their column densities can be derived without knowledge 
of the temperature or density of the absorbing gas. By using both the o-H^^O and 0-H2O 
absorption lines, we are able to constrain the column- averaged 0-H2O abundances in each 
fine-of-sight cloud to within about an order of magnitude. We find N(o-H20)/N(H2) = 
8.1 X 10~^ — 4 X 10~^ for three clouds with optically thin water fines. In three additional 
clouds where the H2O lines are saturated, we have used observations of the Hg^O ground- 
state transition to find upper limits to the water abundance of 8.2 x 10~^ — 1.5 x 10~^. 
The 0-H2O abundances are similar to those determined for a line-of-sight water absorption 
feature towards W51 and Sgr B2 (Neufeld et al. 2000; 2002) but are higher than those seen 
from water emission lines in molecular clouds (Snell et al. 2000b). We measure the OH 
abundance from the average of the 1665 and 1667 MHz observations and find N(0H)/N(H2) 
= 1.2 X 10"^ - 1.1 X 10-6. 

If dissociative recombination is the primary formation pathway for OH and H2O 
then the abundances of these 2 species depends sensitively on the branching ratios ( 
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/h20 ■ /oh '■ fo)- Observations and theoretical work has set these branching ratios to 
fairly discrepant values of either 0.25:0.75:0 or 0.05:0.65:0.3. However, based on a simple 
analytical model, none of the features seen in W49 appear to be consistent with either 
value of the branching ratio. This suggests that a simple model which incorporates only 
dissociative recombination is inadequate, and that additional chemical effects need to be 
considered. Building on the work of Neufeld et al. (2002), our photo-chemistry models 
provide an additional explanation for the observed OH and H2O column densities by 
including depth-dependent photodissociation, neutral-neutral reactions which preferentially 
destroy the OH. The photo-chemistry models can explain the observed OH and H2O column 
densities if /i/jO = 0.25 but not if /h^o — 0.05. These gas-phase models are also able 
to produce results that are consistent with the observed H2O/CO column density ratios, 
without having to resort to freeze-out of CO molecules onto dust grains. However, it is 
possible that atomic oxygen can stick to the grains and rapidly converted to water ice which 
is photodcsorbcd. One attractive feature of this model is that we can still produce high 
H20-to-OH column density ratios in lower density gas (n < 10^ cm~^) with pressures that 
more closely match those observed in other diffuse clouds. 
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Table 1: Spectral Line Observations of W49A 



Telescope 


Line 


On Source Int. Time 


la rms 






(hours) 


(mK) 



SWAS [C 1] -^3 5 

13C0 J = 5^ 4 63 8 

H2O lio loi 63 8 

H^^O lio ^ loi 342 5 

FCRAO i^CO J=1^ 3 50 

13^0 j = 1 ^ 3.2 14 

C^^O J=1^0 4 13 

KOSMA J = 3^2 0.5 60 

Arecibo OH 1665 MHz 0.17 517 

OH 1667 MHz 0.17 673 



Table 2: Gaussian Fit Parameters for the Line-of-Sight Absorption Features:H20 & Hg^O 



Line 


ATX ^ 


Vlsr 


AVfwhm 


JATXdV 




To 




(K) 


(km s~^) 


(km s~^) 


(K km s-i) 






H2O lio - 


> loi -0.067 


33.5 


3.5 


-0.25 


0.22 


1.9 




-0.080 


39.5 


4.1 


-0.35 


0.06 


> 2.7 




-0.051 


53.5 


6.0 


-0.33 


0.40 


1.0 




-0.085 


59.6 


4.0 


-0.37 


0.006 


> 5.1 




-0.073 


63.3 


2.7 


-0.21 


0.15 


>1.9 




-0.032 


68.0 


5.0 


-0.17 


0.63 


0.5 


Hi^O lio - 


loi < -0.012 


34.0 


4.6 


< -0.059 


> 0.86 


< 0.15 




< -0.006 


39.8 


3.2 


< -0.020 


> 0.93 


< 0.07 




< -0.003 


54.0 


5.0 


< -0.015 


> 0.97 


< 0.04 




< -0.007 


59.5 


3.4 


< -0.025 


> 0.92 


< 0.09 




< -0.003 


63.3 


2.7 


< -0.010 


> 0.97 


< 0.04 




< -0.003 


69.0 


5.1 


< -0.015 


> 0.97 


< 0.04 



1 — ATX denotes the baseline subtracted antenna temperature which, for absorption lines, 
is a negative quantity. 
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Table 3: Gaussian Fit Parameters for the Line-of-Sight Emission Features: [C I] 



Line 


ATI 


Vlsr 


AVfwhm 


JATXdV 


/AT^rfl^ (corrected)^ 




(K) 


(km s"""^) 


(km s~^) 


(K km s"^) 


(K km s-i) 


[C I] - 


Po 0.11 


34.0 


4.6 


0.56 


1.00 




0.45 


39.8 


3.2 


1.52 


1.26 




0.13 


54.0 


5.0 


0.71 


0.65 




0.36 


59.5 


3.4 


1.35 


0.66 




0.62 


63.3 


2.7 


1.76 


2.90 




0.12 


69.0 


5.1 


0.70 


0.44 



^ - Corrected to a a 1.5' beam (see Section 4.2). 
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Table 4: Gaussian Fit Parameters for the Emission Features: CO Isotopologues 



Line 








ATI 
(K) 


^LSR 

(km s~^) 


(km s~^) 


/ ATXdV 
(K km s-i) 




J = 


1 


^ " 


0.41 


33.5 


2.0 


0.88 










1.43 


39.3 


3.0 


4.57 










0.33 


54.0 


3.7 


1.33 










0.76 


59.0 


3.3 


2.62 










1.47 


63.0 


2.9 


4.64 










0.36 


69.0 


2.4 


0.94 




J = 


3 


^ 2 


1.42 


39.3 


3.0 


4.61 










0.24 


54.0 


3.5 


0.88 










0.59 


58.5 


3.3 


2.08 










1.49 


63.0 


2.7 


4.22 










0.31 


69.0 


5.1 


1.68 




J = 


1 


^ " 


0.03 


33.5 


2.5 


0.09 










0.68 


39.4 


1.7 


1.23 










0.02 


53.5 


6.0 


0.11 










0.13 


59.4 


3.2 


0.45 










0.96 


63.1 


1.5 


1.52 










0.03 


68.7 


4.6 


0.15 




J = 


5 


^4 


<0.008 










J = 


1 


^ " 


0.09 


39.4 


1.0 


0.09 










0.18 


63.1 


1.0 


0.19 



" - data smoothed to 1.5' resolution. 
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Table 5: Gaussian Fit Parameters for the Line-of-Sight Absorption Features: OH 



Line 


AT4 


Vlsr 


AVfwhm 


/ ATXdV 


Fi/Fc 


To 




(K) 


(km s ) 


[km s ) 


(K km s"^) 






OH 1665 MHz 


-1.3 


33.6 


3.5 


-4.8 


0.995 


0.005 




-0.0 


oy.o 


1 7 
1. / 


-10.41 


n Qfi7 






-1.9 


54.2 


4.7 


-9.7 


0.992 


0.008 




-3.9 


59.8 


3.5 


-14.5 


0.984 


0.016 




-14.2 


62.9 


1.8 


-27.1 


0.943 


0.059 




-1.0 


68.8 


2.8 


-3.0 


0.996 


0.004 


OH 1667 MHz 


-3.9 


33.6 


3.4 


-14.3 


0.985 


0.016 




-8.1 


39.3 


2.4 


-20.8 


0.968 


0.033 




-3.9 


54.2 


4.7 


-19.6 


0.985 


0.016 




-12.0 


59.8 


2.9 


-37.6 


0.952 


0.049 




-24.6 


62.9 


1.8 


-46.4 


0.902 


0.103 




-1.5 


(38.8 


2.8 


-1.5 


0.991 


0.006 



Table 6: LVG Model Results 





log(nj/2) 


log(N(i3CO)) 




log(N(C°)) 


N(C<')/N(i2co)» 


(km s~"^) 


(cm~^) 


(cm~^) 


(mag) 


(cm~^) 




33.5 


3.48 


14.16 


0.1 


16.48 


3.8 


39.4 


3.52 


15.44 


1.9 


16.99 


0.6 


53.5 


3.48 


14.37 


0.2 


16.59 


3.0 


59.4 


3.48 


14.92 


0.6 


16.91 


1.8 


63.1 


3.18 


15.70 


3.4 


17.20 


0.6 


68.7 


3.48 


14.45 


0.2 


16.56 


2.3 



a Assuming N(i2CO):N(i3CO) = 55:1. 

b The visual extinction in the CO emitting layer. Add approximately 2 magnitudes to get 
the total Ay. 
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Table 7: Estimated Water Column Densities and Abundances in W49A 









Vlsr (km s ^) 










33.5 


39.5 


53.5 


59.6 


63.3 


68.0 


To(o - H^^O)AV (km s"^) 


< 0.69 


< 0.23 


< 0.18 


< 0.29 


< 0.10 


< 0.18 


N{o-H^^O) (cm-2) 


< 3.2 X IOI2 


< 1.1 X IOI2 


< 8.4 X 10" 


< 1.3 X 10^2 


< 4.5 X 10" 


< 8.6 X 10" 


N(o-H20) " (cm-2) 


< 1.6 X IQis 


< 5.4 X 10^* 


< 4.2 X 10^* 


< 6.7 X lO^" 


< 2.3 X lO^'* 


< 4.3 X 10^" 


To{o- H20)AV (km s-i) 


6.5 


>11.3 


6.2 


> 20.6 


> 5.2 


2.6 


N{o- H2O) (cm-2) 


3.2 X IQi^ 


> 5.5 X 10^^ 


3.0 X 10^3 


> 1.0 X lO^"* 


> 2.5 X 10^3 


1.3 X 10" 


N{CO) (cm-2) 


8.0 X lO^s 


1.5 X 10^'^ 


1.3 X 10^8 


4.6 X lOi'* 


2.8 X 10^'^ 


1.6 X 10^^ 


Af(() - IhOlNllL^) 


< 2.0 X 10-'' 


< 3.6 X 10-"^ 


< 3.3 X 10-'' 


< 1.5 X 10-« 


< 8.2 X 10-** 


< 2.8 X 10-*' 


iV(<) - Ih())/y(ll2) '' 


1.0 X 10 " 


> :-!.() X 10 


2.3 X 10 ~ 


> 2.2 X 10 ' 


^ 9.2 X 10 


8.1 X 10 



a Assumes N (Hl^ O) / N (Hf O) = 500. 

& Prom LVG calculations (this paper) assuming N(i2cO):N(l3CO) = 55:1. 
c Prom the Hj^O results assuming that CO:H2 = lO"*. 
d Prom the H2O results assuming that CO:H2 = 10"*. 



Table 8: Estimated OH Column Densities and Abundances in W49 





Vlsr (km s ^) 






54.2 


68.8 


ToAV{1665MHz) (km s'^) 


0.04 


0.01 


ToAV{1667MHz) (km s"^) 


0.08 


0.02 


N{OH) 


1.4 X 10^^ 


3.5 X 10^3 


N{OH)/N{H2) 


1.1 X 10-6 


2.3 X 10-^ 


N{H,0)1JN{0H) 


0.26 


0.43 



" - Prom Table 7, row 5 x 4/3 to account for the ortho/para ratio. 
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Fig. 1. — 450 fjLm continuum emission from W49A from SHARC at the CSO with a beamsize 
of approximately 9". The peak flux is ~ 700 Jy but the image scale is truncated at 300 Jy 
to better show the weak, extended emission. 




Fig. 2. — Comparison of the 3 spectral lines detected by SWAS in W49A. (top) Iiq-Ioi- 
(middle) [C I] =^Pi-^Po. (bottom) ^^CO J = 5 ^ 4. 



W49 




Fig. 3. — Comparison of the FCRAO and KOSMA observations of W49A. All spectra have 
been smoothed to the angular resolution of the SWAS ^^CO J = 5 — > 4 hne (~ 3.8'). 
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Fig. 4. — Line to continuum flux ratio in W49A for the SWAS observations of (top) H2^0 
lio ~ loi and (bottom) H2^0 lio — loi- To account for sliglit tilts and curvatures in the 
raw spectra, 1st order baselines were subtracted from the H2^0 and H^^O spectra, prior to 
calculating the line to continuum flux ratios. 
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Fig. 5. — Line to continuum flux ratio in W49A for the Arecibo observations of the (top) 
OH 1667 MHz and (bottom) OH 1665 MHz hues. 
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Fig. 6. — Comparison of observations/limits with results of PDR models. Observations are 
indicated by a black dot or line labeled with the OH velocity of the feature. Model results are 
shown for clouds with Go = 1.7 and Ay— I (sohd curve), 2 (dashed curve), 3 (dotted curve), 
and 4 (dash-dotted curve). For each Ay, results are shown for cloud densities of 10^, 10^ and 
lO"' cm^^, from left to right respectively along each curve. A single point (labeled as "grain") 
indicates the computed H2O and OH column densities from a model with n = 10^ cm^^. 
Ay = 3, and the catalytic formations of water on grain surfaces followed by photodesorption. 
Straight dashed lines follow from the analytical model of Neufeld et al. (2002), where the 
upper line is for branching ratio /H2O '■ /oh : /o of 0.05:0.65:0.30 and the lower line is for 
0.25:0.75:0.0. 
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Fig. 7. — Comparison of the H2O and CO observations with results of PDR models. 
Observations are indicated by the horizontal lines indicating the upper and lower limits of 
the H2O column densities as determined from the H2^0 and H2O observations respectively, 
and are labeled with the H2O velocity of the feature. Model results are shown for clouds 
with Go = 1.7 and k.y=\ (solid line), 2 (dashed line), 3 (dotted line), and 4 (dash-dotted 
line). For each Ay, results are shown for cloud densities of 10'^, 10'^ and 10^ cm~^, from 
bottom to top respectively along each curve. 
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Fig. 8. — Comparison of observed [C I] 492 GHz and ^^CO J = 1 — > emission with results 
of constant density PDR models. Curves show the model results, with the same hne types 
as in the previous figure. Cloud densities are 10^, 10^ and 10^ cm~^ and increase from left 
to right. Points indicate the observed [C I] and ^^CO integrated intensites, corrected as in 
the text, for the line-of-sight toward the 450//m continuum source (with a 1.5' source size). 



